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ABSTRACT 

We present high-resolution spectroscopic observations of the Li I resonance 
line in a sample of 62 stars that belong to 31 common-proper motion pairs with 
twin F or G-type components. Photospheric abundances of lithium were derived 
by spectral synthesis analysis. For seven of the pairs, we have measured large 
lithium abundance differences. Eleven other pairs have components with similar 
lithium abundances. We cannot determine if the remaining 13 pairs have lithium 
differences because we did not detect the Li I lines, and hence we can only provide 
upper limits to the abundances of both stars. Our results demonstrate that twin 
stars do not always share the same lithium abundances. Lithium depletion in 
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solar-type stars does not only depend on age, mass, and metallicity. This result 
is consistent with the spread in lithium abundances among solar-type stars in 
the solar-age open cluster M67. Our stars are brighter than the M67 members 
of similar spectral type, making them good targets for detailed follow-up studies 
that could shed light on the elusive mechanism responsible for the depletion of 
lithium during the main-sequence evolution of the Sun and solar-type stars. 

Subject headings: stars: abundances — stars: late-type - binaries: visual - stars: 
evolution - stars: fundamental parameters 

1. Introduction 

Lithium has a rich life. It is created in many environments (primordial nucleosynthe- 
sis, cosmic ray spallation in the interstellar medium, spallation and fusion reactions around 
compact relativistic objects, flares, and red giant thermal pulses), and it is destroyed in the 
interior of stars by collision with protons at a temperature of about 2.7xl0 6 K. The Sun 
has a lithium abundance more than two orders of magnitude lower than meteoritic material 
(Miiller, Peytremann, & de la Reza 1975), which reflects the composition of the presolar neb- 
ula. The Sun probably did not deplete lithium during the pre main-sequence evolution, but 
rather during the slow main-sequence evolution, although the exact mechanism has not been 
identified yet (Martin 1997, 1998). Recent models of lithium depletion in solar-type stars 
suggest that two types of mixing may be at work during main-sequence evolution; namely 
overshooting and rotational mixing (Umezu & Saio 2000). The history of rotational mixing 
may be strongly influenced by the initial conditions of angular momentum distribution in 
the protoplanetary disk. Thus, it is conceivable that there could be a connection between 
the rate of lithium depletion by rotational mixing and the presence of companions (stellar 
or substellar) to the stars. 

Binary stars with twin components are interesting for understanding binary formation 
and evolution. Common proper motion (CPM) pairs of twins allow us a test of the validity of 
stellar evolution models because their separations are so large (>100 AU) that each member 
has probably evolved independently. If all the properties of stars are determined by age, 
mass, and chemical composition, twins should be identical. If the lithium abundance of 
a main-sequence star is determined by its age, mass, and chemical composition, pairs of 
twins should show the same lithium abundances. In the early work of Herbig (1965), it 
was already apparent that there was a scatter in lithium abundances among pairs of twins. 
Herbig's study included 53 UMa A and B, where lithium was detected in A but not in B; 
16 Cyg A and B, where lithium was detected in A but not in B; and 53 Aqr A and B, where 
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lithium was detected in both stars. 

The system 16 Cyg has received some attention recently. The two stars of the pair are 
nearly identical. The primary has V=5.96, T cff =5785 K and log g=4.28. The secondary has 
V=6.20, T cff =5747 K and log g=4.35. Despite those similarities, the members of this pair 
have two fundamental differences. King et al. (1997) have shown that the primary has a 
lithium abundance higher than that of the Sun (logN(Li) = 1.27±0.05 in the customary scale 
of log N(H)=12), while the secondary has a much lower lithium abundance (logN(Li)<0.60). 
The large difference in lithium abundance is surprising because both stars should be the same 
age and have very similar masses. Another significant difference between the two stars is that 
the secondary harbors a giant planet with a minimum mass of 1.5 Jupiters, an eccentricity 
of e=0.63 and a period of 800.8 days (Cochran et al. 1997). On the other hand, no giant 
planet has yet been detected around 16 Cyg A. However, such nondetection does not rule 
out the presence of giant planets with long periods or smaller planets with short periods. 

Gonzalez & Laws (2000) and Ryan (2000) have studied the distribution of lithium 
abundaces among planet-harboring stars, and have reached opposite conclusions. While 
Gonzalez & Laws claim that lithium is overdepleted in stars with giant planets, Ryan claims 
that it is normal. The disagreement comes from the choice of stars to be compared with the 
host stars of extrasolar planets. This is a difficult task because the ages of most field stars 
are not well known. 

Pairs of twins offer the possibility of comparing stars that are coeval, and have nearly 
the same mass and metallicity. Thus, we may expect that the systematic study of binaries 
may shed light on whether there is any connection between the presence of giant planets 
and lithium depletion. Gratton et al. (2001) have recently reported an abundance analysis 
for six wide binaries. They find that the primary star of the binary HD 219542 has higher 
iron and lithium abundances than the secondary star, which leads them to suggest that the 
primary has ingested a planet. 

In this paper we report lithium abundances for a sample of 31 visual binaries. The tar- 
gets have been chosen to be as similar as possible to 16 Cyg in spectral type and magnitude. 
The paper is organized as follows: In Section 2 we describe the observations. In Section 3 we 
present the abundance analysis. In Section 4 we discuss each pair of twins individually when 
we could detect the Li I line in at least one of the stars, and we compare them with lithium 
abundances in open clusters and theoretical models. Section 5 contains the discussion, and 
in Section 6 we present our conclusions. 
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2. Observations 

We selected our targets from the list of CPM pairs compiled by Halbwachs (1986). We 
filtered the sample by requiring that the difference in V magnitudes should be less than 0.4, 
and the spectral type should be between F8 and K0. The declination of the stars had to be 
north of —30 degrees. Our sample is listed in Table 1. 

Observations were obtained during several observing runs at Lick and Keck observa- 
tories between January 1998 and March 2001. The binaries BD+08 4386, BD+69 0993, 
BD+74 0718, and BD+76 0835 were observed at Keck. All the other binaries were observed 
at Lick. The observing log is given in Table 1. 

The Lick observations were carried out using the Hamilton echelle spectrograph on the 
Shane 3-meter telescope. This instrument provides 92 orders covering the spectral range 
from 490 nm to 890 nm at a resolution of 48,000. A standard reduction was performed 
(Valenti 1994), which included background subtraction, cosmic ray removal, and flat fielding. 
Wavelength calibration was made using ThAr lamp exposures. 

The Keck observations were obtained with the HIRES echelle spectrograph (Vogt et al. 
1994). The resolution was 67,000 and the coverage was from 600 to 850 nm. Data reduction 
was performed using IDL routines similar to those used for the Lick data. 

3. Analysis 

We have compiled the available photometry in the literature for our program stars and 
we provide it in Table 2. The V magnitudes, B — V colors, and spectral types were taken from 
SIMBAD, except for the particular cases listed in Section IV, and for the notes to Table 2. 
This photometric database is rather heterogeneous. Some of our stars are included in the 
2MASS second incremental release point source catalog (PSC). Their 2MASS J — K colors 
are included in Table 2. Temperatures were derived using the empirical scale of Alonso, 
Arribas & Martinez- Roger (1996) for solar metallicity. The temperatures obtained from the 
J — K colors are given in Table 3. Typical error bars in J — K colors are ±0.03, corresponding 
to ±180 K in T e fj. These temperatures were derived to compare with our, presumably more 
reliable, temperatures obtained from synthetic spectra fitting to the observed spectra in the 
lithium region. This comparison is shown in Figure 1 for 26 stars that have available 2MASS 
JK photometry. The mean difference between the T efr derived from J — K and from fits 
is —65 K and the standard deviation is 270 K. We conclude that there is no significant 
systematic deviation of the T e g derived from J — K colors with respect to those derived from 
synthetic fits. 
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For each spectrum we determined the continuum level using the DECH20 software 
(Galasutdinov 1992). To detect the real continuum level, we smoothed the observed spectra 
using three loops of a 3-point smoothing boxcar. Then the continuum points were fitted 
by splines to use in the fine analysis procedure. The S/N ratios per pixel given in Table 1 
were obtained by measuring the standard deviation of the normalized continuum in the 
spectral range 669.3-669.5 nm. We carried out the computations using the LTE spectral 
synthesis program WITA6 (Pavlenko 2000), which is a modified version of the program 
used by Pavlenko et al. (1995). Our computations of synthetic spectra were carried out 
using plane-parallel model atmospheres in LTE, with no energy divergence. The model 
atmospheres were taken from Kurucz's grid available on the web. 

We used metallicities in the range [m/H]=0.0 to [m/H]=-0.2 for fitting the observed 
spectra 1 except for HD 54046 and HD 54100, which are known metal poor stars (Norris 
1986). Chemical equilibrium was computed for the mix of ^100 molecular species, with 
data taken from Tsuji (1973). ATLAS9 (Kurucz 1993) grid of continuum opacity sources 
were used in our computations. We used VALD (Piskunov et al. 1995) for the atomic line 
list and damping constants. For some lines the damping data are missing; in that case we 
used Unsold (1955) approach to compute them. We had to increase the log gf of the Fel 
line at 670.510 nm by a factor of 2 (from 0.0319 given by VALD to 0.0640) to find good 
fits to the observed spectra of the Sun as a star (Kurucz et al. 1984). Computations were 
carried out for a fixed microturbulent velocity of V t —2 km s _1 and a Voigt profile for every 
absorption line. 

Synthetic spectra were computed with a step of 0.002 nm. Instrumental profile and 
macroturbulence broadening effects were modeled by convolution with a Gaussian of half- 
width in the range 0.010 nm to 0.018 nm. For most stars we could desribe the profile 
with instrumental broadening only, without any rotation. That provides upper limit of 
vsim=3 km s -1 . No enhanced van-der Waals broadening was used. For two stars rotational 
broadening was taken into account (HD 6872 A and BD+65 1044; vsim=14 km s _1 ). For 
each star we found the best-fitting synthetic spectrum using a least-squares minimization 
algorithm. Figures 2, 3 and 4 show typical examples of synthetic fits to the echelle spectra. 

For each LTE lithium abundance, we determined NLTE abundance corrections usind 
the grid of curves of growth of Pavlenko & Magazzu (1996). In Table 3 we summarize our 
results. We give the temperatures, metallicities and lithium abundances (including NLTE 
corrections) that we used in the best fits to the observed spectra. Equivalent widths of the Li I 
line measured by direct integration of the line profile in the observed spectrum are given in 



1 We refer to the solar abundances given by Grevesse (1984) using the convention [m/H]=0 
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Table 3 as well. They can be used to derive lithium abundances from curve of growth (COG) 
computations. For the stars with 2MASS photometry we have obtained COG Li abundances 
using the temperatures derived from the (J — K) — T e g- calibration of Alonso et al. (1996) and 
the WITA code. Those abundances are listed in Table 3 for comparison with the abundances 
obtained from the synthetic spectra. Figure 5 illustrates this comparison for the same 26 
stars that were used in Figure 1. The mean difference between the lithium abundances 
derived using these two methods is log N (Li) cog - log N (Li) fit = —0.08 dex, and the standard 
deviation is 0.17 dex. We conclude that there is no significant systematic discrepancy between 
the two methods of obtaining lithium abundances. We adopt the uncertainty of 1 cx=0.17 dex 
in our absolute lithium abundances. In our computations we have assumed that our stars 
are dwarfs. A variation of log g ±0.5 dex provides a maximum AlogN(Li)<0.1 dex. 



4. Comments on Individual Stars 

The main goal of this paper is to determine whether stars with the same age and initial 
chemical composition, and nearly equal masses, have depleted lithium at the same rate. We 
limit our discussion to pair of twins for which we could detect the Li I line in at least one of 
the stars. When we could not detect the Li I line in any of the stars, we could not tell if the 
lithium abundances are different or similar, and therefore those pairs are not discussed any 
further. After a literature search using SIMBAD, we found that some of the pairs deserve 
individual comments: 

• HD 6872 A and B: Oblak & Chareton (1980) and Duncan (1984) included HD 6872 A 
and B among a sample of stars that are sufficiently evolved for the ages to be estimated 
from T eff and absolute magnitudes using the displacement from the zero-age main- 
sequence (ZAMS). Both groups obtained uvby/3 photometry and derived T eff of 6380 K 
(Oblak & Chareton) and 6300 K (Duncan) for the primary. For the secondary, they 
gave 6251 K (Oblak & Chareton) and 6084 K (Duncan). The discrepancy in T e g 
is mainly due to their use of different equations of b — y index versus temperature. 
Soderblom, Duncan & Johnson (1991) converted the b — y color of Duncan to B — V 
color. In Table 1, we adopt the Soderblom et al. color, which is different from the 
color provided by SIMBAD. Our best fit temperatures are 6250 K for both stars. We 
find good synthetic fits to our spectra with the same T eff for both stars (Figure 3). 
The fainter star has a higher lithium abundance than the hotter star by at least 0.8 
dex. Adopting the T cff of either Oblak & Chareton or Duncan does not change the 
conclusion that the lithium abundances of these two stars are significantly different. 
The primary star has broader lines than the secondary. Our best fit to the spectrum 
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of HD 6872 A has rotational broadening of vsim=14 km s 1 . 

• HD 39274 and HD 39275: The lithium abundance difference between these two stars 
is O.f dex, which is less than the 1 a uncertainty in our analysis. HD 39274 is a binary 
with separation 0.225 arcsec (McAlister et al. 1989). 

• HD 54046 and HD 54100: Metal poor stars. Norris (1986) reported [Fe/H]=-0.6 and 
Eggen (1998) gave [Fe/H]=-0.45 based on photometric colors. 

• HD 8624 and HD 8610: Tokovinin (1999) found that HD 8624 is a double-lined spec- 
troscopic binary with a period of 14.9 days. HD 8610 has a K-type companion, and 
therefore the system is a quadruple. Tokovinin gave V=8.83 for both components of 
the HD 8624 close system, and inferred a spectral type of G3. 

• HD 98744 and HD 98745: HD 98745 has a lithium abundance 1 dex or more higher 
than HD 98744. Duncan (1984) found T eff =6050 K and 6039 K for HD 98744 and HD 
98745, respectively. Our best fit T c g is 6250 K for both stars. Adopting the T c g of 
Duncan does not change the conclusion that the lithium abundances of these two stars 
are significantly different. 

• HD 167215 and HD 167216: In Table 2 we give the B - V colors given by SIMBAD, 
because they are more consistent with the temperatures that we find in our analysis, 
than the B — V colors of Soderblom et al. (1991). These latter authors derive T e g of 
5984 K and 6040 K for HD 167215 and HD 167216, respectively. We find that we need 
hotter temperatures to fit the spectra in the lithium region. We also need subsolar 
metallicities. The lithium abundance difference is 0.1 dex. The Li I resonance line 
of HD 167215 is clearly stronger than that of HD 167216. According to Soderblom 
et al, HD 167215 is more evolved above the ZAMS than HD 167216. The gravity of 
HD 167215 may be somewhat lower than that of HD 167216, and this could explain 
the difference in Li I line strengths. 

• HD 224984 and HD 224994: These two stars have nearly identical spectra. HD 224994 
itself is a close binary (separation=0.24 arcsec in 1988; McAlister, Hartkopf, & Franz 
1990), so the system is at least triple. The two stars have very similar spectra, including 
the Li I lines. We obtain lithium abundances of logN(Li)=2.3, which is about 0.8 dex 
below the cosmic value. Therefore, this is an example of a pair of twins that have 
depleted the same amount of lithium. The presence of a low-mass star companion at 
about 20 AU of HD 224994 has not changed the lithium depletion rate of this star with 
respect to HD 224984. 
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5. Discussion 

Lithium abundance is very sensitive to the mass (and thus T e g) of stars. In the temper- 
ature range of our program stars (6500-5250 K), lithium abundances generally decrease with 
decreasing T e g, although at the hot end of this temperature range, the F-type lithium dip 
starts to kick in (Boesgaard & Tripicco 1987). If the lithium abundances of the components 
of pairs of twins were dominated by the mass dependance of lithium depletion, we would 
expect a correlation between the difference in T cff and the difference in logN(Li), in the sense 
that the cooler components should show systematically lower lithium abundances than the 
hotter ones. Figure 6 shows that such correlation is not present in the data. We conclude 
that the observed pattern of lithium abundances in pairs of twins cannot be explained solely 
with the dependance of lithium depletion on stellar mass. 

We have found five pairs with nearly identical T e ff and nearly identical logN(Li), and 
four pairs with nearly identical T e ff and very different logN(Li). Hereafter we call these 
latter four pairs "16 Cyg analogs", namely: HD 6872 A/B, BD+60 0269/BD+60 0271, 
BD+13 2311 A/B, and HD 98744/HD 98745. The common characteristics of" 16 Cyg analogs" 
is that they are pairs of twins with lithium abundance difference between the two stars of the 
pair that exceeds 0.5 dex (which corresponds to 3 a significance in our synthetic fit analysis) 
and difference in T cff <200 K. 

Figure 7 illustrates the comparison between the pairs of twins and the Hyades cluster in 
a T c fj versus logN(Li) diagram. The Hyades single star members define a tight relationship 
between T eff and logN(Li). In general the twins do not follow the Hyades relation. A possible 
explanation for this discrepancy could be that there are many tidally locked binaries (TLB) 
hidden among the twins. We do not know of any example of a TLB in our sample. 

HD 8610 has a lithium abundance similar to Hyades members with the same T e ff. The 
lithium- T eff locus in the Hyades declines steeply from 5500 K to 5250 K, so one would expect 
that HD 8624 should have a lithium abundance much lower than that of HD 8610, which 
is contrary to what we have found. HD 8624 is a double-lined spectroscopic binary with a 
period of 14.91 days (Tokovinin 1999) and an eccentricity e=0132. Even though it is not 
strictly a TLB because the orbit is not completely circularized, it may have preserved lithium 
in a manner similar to the TLBs studied by Barrado y Navascues et al. (1997), which show 
a trend of higher lithium abundances when compared to single stars. However, TLBs are 
relatively rare among the general population of solar-type stars. A program to monitor the 
radial velocity of 10 of our pairs has been started using HIRES with the iodine cell. There 
are already enough data for 16 of the stars to find TLBs, but none has been seen (Marcy 
2002, private communication). We conclude that it is very unlikely that TLBs can account 
for most of the discrepancy between the pairs of twins and the Hyades. 
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The average age of the pairs of twins is probably larger than the Hyades cluster age 
(~600 Myr). It is informative to compare the twins with the well-studied solar-age cluster 
M67. In Figure 7 we compare the pairs of twins with the members of M67. The M67 stars 
show a large spread of lithium abundances at a given T e g, which is similar to the behavior 
exhibited by our binaries. The evolution of lithium from the Hyades to M67 implies a 
mechanism of lithium depletion during the main-sequence lifetime of solar-type stars that 
is not very mass dependent. We have not found in the literature any models that can 
explain the evolutionary pattern of lithium abundances implied by the comparison between 
the Hyades and the M67 solar-type stars. We conclude that there is currently not a good 
theoretical explanation for the lithium abundances observed in our pairs of twins and the 
open clusters. In future papers we plan to tackle this problem observationally via searches 
for connections between the lithium abundances in the pairs of twins and the presence of 
companions (stellar or substellar), chromospheric activity, rotation, and the abundances of 
other chemical elements. At the moment it would be premature to suppose that the presence 
of planets had anything to do with the spread in lithium abundances. 



6. Final Remarks 

We report high-resolution spectroscopy observations of 62 solar-type stars in 31 common 
proper motion pairs of twins. These binaries are sufficiently separated that the evolution 
of each star is not expected to be influenced by its stellar companion. Lithium abundances 
have been derived from theoretical fits to the observed spectra using model atmospheres. 
A comparison with COG abundances derived from Li I equivalent widths and temperatures 
obtained from near-infrared colors shows that the absolute lithium abundances are reliable 
to within 0.17 dex (1 a). 

We have found seven pairs of twins with large lithium abundance differences. The pat- 
tern of lithium abundances in our sample of twins seems to be midway between those of the 
Hyades (tight lithium-T efr relation) and M67 (large lithium spread for any T eff ) open clus- 
ters. We have not found any evolutionary models in the literature that provide a satisfactory 
explanation for the observed spread in lithium abundances in solar-type stars. We plan to 
follow up with additional studies of our sample of twins to obtain an overall understanding 
of all the effects that may come into play in the evolution of solar-type stars. 

This paper is based on data collected with the 3-m Shane telescope at Lick Observatory, 
California, and the Keck I telescope at the W. M. Keck Observatory, Mauna Kea, Hawaii. 
The Keck observatory is operated as a scientific partnership among the California Institute 
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istration. The Observatory was made possible by the generous financial support of the W. 
M. Keck Foundation. This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. The authors wish to extend special thanks to those of Hawaiian 
ancestry on whose sacred mountain of Mauna Kea we are privileged to be guests. Without 
their generous hospitality, the Keck telescope observations presented therein would not have 
been possible. We thank Louise Good for correcting the English of the manuscript, and 
Geoff Marcy for communicating results of the radial velocity program for a subset of our 
sample. We are also indebted to Ann Boesgaard and David Barrado y Navascues for sending 
data files for lithium abundances of M67 stars, and Nina Kharchenko for astrometric data. 
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Fig. 1. — Comparison of effective temperatures obtained from J — K colors and from the 
best synthetic fits to the observed spectra. 
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Fig. 2. — Solid lines: Lick Hamilton echelle spectra of the CPM pair HD 4552 A and B. 
Dashed lines: Synthetic spectra for T c ff=6000 K, [m/H] =— 0.2, and logN(Li) = 1.6 (above), 
and for T cff =5750 K, [m/H] = -0.2, and logN(Li)=2.3 (below). 
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Fig. 3. — Solid lines: Lick Hamilton echelle spectra of binary HD 6872 A and B. 
Other name for HD 6872 B is BD+12 0090. Dashed lines: Synthetic spectra for 
T cff =6250 K, [m/H]=-0.2, and logN(Li)=1.8 (above), and for T cff =6250 K, [m/H]=-0.2, 
and logN(Li)=2.6 (below). 
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Fig. 4. — Solid lines: Lick Hamilton echelle spectra of the CPM pair HD 9518 A 
(BD+60 0269) and B (BD+60 0271). Dashed lines: Synthetic spectra for T cff =6250 K, 
[m/H]=0.0, and logN(Li)=2.55 (above), and for T cff =6250 K, [m/H]=0.0, and logN(Li) = 1.8 
(below). 
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Fig. 5. — Comparison of lithium abundances obtained with J — K colors and COG calcula- 
tions and with the best synthetic fits to the observed spectra. 
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Fig. 6. — Differences in the T e fr adopted for pairs of twins versus the differences in the lithium 
abundances derived by us (solid symbols) and by previous investigations (open symbols). 
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Fig. 7. — Lithium abundances versus T e g for pairs of twins. Our lithium measurements are 
given as solid dots, with the components of each binary linked by a solid line. The effective 
temperatures of some pairs are shifted by 30 K to avoid confusion. The components of 16 Cyg 
are denoted with an open dot and a solid line. The locus of lithium abundances in single 
stars members of the Hyades cluster is shown with a dotted line. The Hyades data have been 
taken from Thorburn et al. (1993). Lithium abundances for M67 members are denoted with 
open hexagons for detected Li I lines, and asterisks for upper limits. Data for M67 stars come 
from Hobbs & Pilachowski (1986), Spite et al. (1987), Garcia Lopez, Rebolo, & Beckman 
(1988), Deliyannis et al. (1994), Pasquini, Randich, & Pallavicini (1997), Jones, Fischer, & 
Soderblom (1999), Barrett et al. (2001), and Barrado y Navascues, & Balachandran (2002). 
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Table 1. Observing Log. 



HD/BD 


Telescope 


Date 


Texp (s) 


S/N 


HD 4552 


3-m Lick 


1998 


Nov 23 


1500 


140 


BD+12 0090 


3-m Lick 


1998 


Nov 23 


1800 


105 


HD 6872 A 


3-m Lick 


1998 


Sept 4 


1000 


140 


HD 6872 B 


3-m Lick 


1998 


Sept 4 


1100 


130 


HD 8624 


3-m Lick 


1998 


Sept 3 


600 


105 


HD 8610 


3-m Lick 


1998 


Sept 3 


600 


130 


BD+60 269 


3-m Lick 


1998 


Sept 3 


800 


110 


BD+60 0271 


3-m Lick 


1998 


Sept 3 


1400 


100 


HD 31208 


3-m Lick 


2001 


March 8 


1800 


150 


BD+07 754 


3-m Lick 


2001 


March 8 


1800 


190 


HD 39274 


3-m Lick 


2001 


March 8 


2500 


190 


HD 39275 


3-m Lick 


2001 


March 8 


2200 


175 


HD 54046 


3-m Lick 


1998 


Jan 22 


1800 


280 


HD 54100 


3-m Lick 


1998 


Jan 22 


1800 


340 


BD+28 1698 


3-m Lick 


1998 


Jan 21 


1200 


210 


BD+28 1697 


3-m Lick 


1998 


Jan 21 


1200 


200 


BD+15 2080 


3-m Lick 


1998 


Jan 22 


2400 


210 


BD+15 2079 


3-m Lick 


1998 


Jan 22 


1800 


250 


BD+34 2091 A 


3-m Lick 


2001 


March 9 


3000 


90 


BD+34 2091 B 


3-m Lick 


2001 


March 9 


3600 


95 


HD 92222 A 


3-m Lick 


2001 


March 8 


3600 


140 


HD 92222 B 


3-m Lick 


2001 


March 8 


3600 


135 


BD+13 2311 A 


3-m Lick 


2001 


March 9 


5130 


225 


BD+13 2311 B 


3-m Lick 


2001 


Jan 22 


4400 


220 


HD 98744 


3-m Lick 


2001 


March 8 


4000 


175 


HD 98745 


3-m Lick 


2001 


March 8 


4500 


175 


HD 111484 A 


3-m Lick 


1998 


Jan 21 


1800 


200 


HD 111484 B 


3-m Lick 


1998 


Jan 21 


1800 


280 


HD 124257 A 


3-m Lick 


1999 


May 26 


2400 


180 


HD 124257 B 


3-m Lick 


1999 


May 26 


1800 


150 


HD 124913 A 


3-m Lick 


1999 


May 27 


1800 


110 


HD 124913 B 


3-m Lick 


1999 


May 27 


1800 


105 


BD+02 2820 A 


3-m Lick 


1999 


May 26 


2000 


95 


BD+02 2820 B 


3-m Lick 


1999 


May 26 


2400 


120 


BD+35 2576 A 


3-m Lick 


1999 


May 27 


2200 


135 


BD+35 2576 B 


3-m Lick 


1999 


May 27 


2000 


100 


BD+13 2830 A 


3-m Lick 


2001 


March 8 


2000 


65 


BD+13 2830 B 


3-m Lick 


2001 


March 8 


2400 


60 


BD+15 2867 B 


3-m Lick 


1999 


May 28 


2200 


80 


BD+15 2867 C 


3-m Lick 


1999 


May 28 


2400 


85 


HD 155674 


3-m Lick 


1999 


May 28 


2000 


85 


BD+54 1862 


3-m Lick 


1999 


May 28 


2000 


110 


BD+65 1043 


3-m Lick 


1998 


Sept 4 


1800 


105 


BD+65 1044 


3-m Lick 


1998 


Sept 4 


2000 


100 


BD+74 0718 


10-m Keck I 


1998 


Sept 20 


200 


95 


BD+74 0719 


10-m Keck I 


1998 


Sept 20 


200 


150 


HD 167215 


3-m Lick 


1998 


Sept 3 


800 


120 


HD 167216 


3-m Lick 


1998 


Sept 3 


800 


110 


BD+69 0993 A 


10-m Keck I 


1998 


Sept 20 


120 


155 
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Table 1 — Continued 



HD/BD 




Telescope 


Date 




Texp (s) 


S/N 


BD+69 0993 


B 


10-m Keck I 


1998 Sept 


20 


120 


115 


BD+08 4386 


A 


10-m Keck I 


1998 Sept 


20 


180 


105 


BD+08 4386 


B 


10-m Keck I 


1998 Sept 


20 


180 


110 


BD+03 4428 


A 


3-m Lick 


1998 Sept 


4 


1500 


120 


BD+03 4428 


B 


3-m Lick 


1998 Sept 


4 


1800 


110 


BD+76 0835 


A 


10-m Keck I 


1998 Sept 


20 


180 


160 


BD+76 0835 


B 


10-m Keck I 


1998 Sept 


20 


180 


150 


BD+01 4575 


A 


3-m Lick 


1998 Sept 


4 


2500 


80 


BD+01 4575 


B 


3-m Lick 


1998 Sept 


4 


2500 


80 


BD+11 5033 


A 


3-m Lick 


1998 Sept 


3 


1200 


100 


BD+11 5033 


B 


3-m Lick 


1998 Sept 


3 


1500 


105 


HD 224984 




3-m Lick 


1998 Sept 


3 


800 


110 


HD 224994 




3-m Lick 


1998 Sept 


3 


1000 


90 
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Table 2. Photometric Data for Program Stars. 



HD/BD 


V 


B - V 


Sp.T. 


J 


J- H 


J - K 


b-y 


ml 


cl 


Sources 


rlLJ 4iD0Z 


fi fi7 


n ^ 

U.OO 


r o 


7 °,A 
I .04 


fl 99 
U. ZZ 


fl 9Q 

u.zy 










idlj-\- iz uuyu 


o on 
y.zu 


n 

U.OO 




n7 
o.U ( 


n Qn 
U.oU 


fl Q/1 










TT n C070 A 
UU Do (Z A 


c nfi 

O.UO 


n a i 

U.41 


r o 


7 1 
( .ID 


n 1 fi 

U. ID 


n 9*5 
U.zo 


n 901 
u.zy l 


n 1 fio 
u. ioy 


n /i9i 

U.4Z1 


A/f 7fi 
1V1 / D 


Un fifi79 R 
LILJ DO I Z JO 


fi K, 1 * 


n fi7 
u.o / 


r o 


7 ^fl 
/ .OU 


fl 9zt 


fl ^9 
u.oz 


fi ^fi 

U.OOO 


fl 1 ^ 
U. 100 


fl ^7fi 
U.O i D 


M 7fi 
1V1 / D 


LLLJ OOZ4: 


fi m 


n fi7 

U.O 1 


Pfi 








fl AOA 
U.4Z4 


fl 9*3/1 
U. Z04 


fl ^^fl 
U.OOU 


PlQzl 


un o«i n 
nu ooiu 


fi m 


n fin 
u.ou 


r o 








U.4UU 


fl 1 Q1 

u. ly i 


U.4UO 


PlQzl 


Rn_Lf;n 9fio 
du-|-ou zoy 


Q t^fi 
O.OO 


n a s 

U.4o 


r o 
















p. Pi i fin 971 


q nfi 
y.uo 


u.4y 


r o 
















Un Q1 oris 
rlU olZUo 


q on 
o.zu 


n fiQ 
U.oo 


IvU 








n ziofi 
u.4yo 


n QOQ 

U.oyo 


n qio 
U.olz 




JolJ-|-Uf (04 


fi "!Q 
O.OO 


n so 
U.oy 


IV u 








U.010 


n /iO/i 

U.4Z4 


n QH7 
U.oU f 


CQO 

oyo 


Un Q007/I 

UU oyz / 4 


fi OA 

o.yo 


n c;i 


LtU 


7 SO 

( .oy 


n qi 
U.ol 


fl Q/1 
U.o4 










un QQ97c; 


fi QQ 

o.yy 


U.o4 


Pfl 


i .Do 


fi ^fi 

U.OD 


fl AA 










un rx/in/lfi 

nu 04U40 


7 79 
I . ( Z 


fi RA 
U.D4 


Pfi 








fl 

u.ooy 


u. too 


fl 9Qfi 

u.zyo 


no 4 

U04 


un K/iinn 

UU D41UU 


7 79 
1 . / z 


U.OO 


r o 








n Qfio 

U.ODZ 


n 1 c^o 
u. ioy 


n q i fi 

U.OlD 


Plfi/1 
Uo4 




fi 97 
O.Z / 


u.oo 


Pfi 
















pn i oo -1 cq7 

DUi^zo ioy i 


fi 9fi 
o.ZO 


u.ou 


Pfi 
















pn i 1 cr 9fiQfl 
JOlJ-|-10 ZUoU 


y.ui 


n As 

U.Do 


















pp. I i c; 9f|7Q 


q 77 
O. I ( 


u.oy 


Pfl 
















pn 1 Q/i 90Q1 A 

DLj-\-o'-t zuyi /\ 


Q 71 


n ^i 

U.Ol 


Pfl 

UU 


O.DU 


fl 9Q 

u.zy 


fl ^fi 
U.OD 










PPl 1 OA 9f|Q1 P< 

DLj-\-o'-t zuyi -D 


Q 7Q 

y. / y 


fl AQ 

u.4y 




o.oD 


fl 9*3 

u.zo 


fl ^9 
U.OZ 










un Q9999 A 


Q AC\ 

y.4u 


n 70 

U. I u 


Pfl 
VjU 


fi 1 7 

0. 1 ( 


fl ^fi 
U.oo 


fl A^ 
U.40 


fl zl^fi 
U.40O 


fl 9fi1 
U. Zol 


fl ^9fi 
u.ozo 


PlfiA 
1J04 


un Q9999 P 


q ^n 
y.ou 


fl 7Q 

u. * y 


Pfl 

uU 


fi 1 A 
0. 14 


fl ^fi 
U.OD 


fl 49 
U.4Z 


fl zl^fi 
U.40D 


fl 9fifl 
U. ZDU 


fl ^7fi 
U.O i D 


no 4 

1J04 


pn i i q oqi i a 
£>J_J-(- lO Zoll rt. 


O Q9 

y.oz 


U.4D 


r o 


fi 7/1 
O. I 4 


n 9/1 

U.Z4 


n 9^ 
U.zO 










pn i i q 9qi i R 
JolJ-|-10 Zoll -D 


fi fi 
O.O 


U.D. 


r o 


fi ^ 
O.OO 


fl 9^ 

u. zo 


fl 9^ 
U.ZO 










Un QQ7/I/L 
riJJ yo i 4^ 


y. io 


fl zlfi 
U.4D 


Pfl 
uU 


fi 99 

o.zz 


fl 9fi 
U.ZD 


fl ^fl 
U.OU 


fl ^fi 
U.OOD 


fl 1 A9 
U. 14Z 


fl ZL1 Q 

u.4iy 


l\/Tl\A7fi 
1V11V1 t O 


un QQ7/I 

H1J yo ( 40 


O "JO 

y.oy 


fl A7 
U.4 ( 


nn 
LtU 


fi /19 
0.4Z 


n 9^ 
U.zO 


n oo 
U.zy 


U.OvjO 


n 1 £9 
U. lOz 


n Q7Q 

U.o i o 


A/Tl\'T7fi 


un iii /i qa A 

HU 111^404 /\ 


fi 77 
O. / f 




Pfl 














PCQQ 


un 1 1 1 A Q.A P 
liU lil4o4 Jo 


o go 
O.OO 




pn 
LtU 














PQnn 


Un 1 9/1 9^7 A 
UU 1Z4Z0 1 r\ 


Q 7 

y . ( 


n Q- 
U.o. 


p c 


fi i n 
o. 1U 


n 99 

U.ZZ 


n 90 
u.zy 










Un 1 9A9^7 R 
riJJ iz^zo/ Jo 


Q fi 
y .o 


1 1 ■ 

1.1. 


p c; 


fi fifi 

O.DO 


fl 9^ 

u. zo 


U.OO 










Un 1 9/1 Ol Q A 
LLLJ lZ4yio r\ 


o n 
y .u 


1 fl 
l.U 


p c; 


fi f^O 

o.oy 


n 9i 

U. Zl 


n 9/1 

U.Z4 










un 1 9/1Q1 Q R 

nu iz4yio Jo 


O 9 

y .z 


"1 1 
1. 1 


p c; 
LtO 


fi 09 

o.yz 


n 17 
U.l / 


n 9n 
U.zU 










Pn_l_ri9 9Q90 A 
JoU-|-UZ ZoZD r\ 


1 n n 

lu.U 


n o 
u.y 


p c; 
















pn i no ocon p 
dL)-\-\JZ zozU Jo 


1 n n 
1U.U 


n n 
U.y 


P 
LtO 
















BD+35 2576 A 


10.20 


0.67 


KO 
















BD+35 2576 B 


10.77 


0.61 


KO 
















BD+13 2830 A 


9.8 


1.1 


KO 


9.28 


0.34 


0.37 










BD+13 2830 B 


9.9 


1.1 


KO 


9.37 


0.33 


0.43 










BD+15 2867 B 


9.7 


1.1 


KO 


9.15 


0.39 


0.46 










BD+15 2867 C 


9.7 


1.1 


K0 


9.16 


0.40 


0.49 










HD 155674 


8.80 


1.20 


KO 
















BD+54 1862 


9.29 


1.24 


K8 
















BD+65 1043 


8.95 


0.99 


G5 


7.19 


0.46 


0.60 










BD+65 1044 


9.67 


0.43 


F8 


8.84 


0.18 


0.21 










BD+74 0718 


11.1: 


0.87 


GO 
















BD+74 0719 


10.09 


0.77 


G5 
















HD 167215 


8.10 


0.48 


F8 








0.357 


0.148 


0.373 


D84 


HD 167216 


8.23 


0.48 


F8 








0.346 


0.160 


0.357 


D84 


BD+69 0993 A 


10.1 


1.0 


G 
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Table 2 — Continued 



HD/BD V B-V Sp.T. J J-H J-K b-y ml cl Sources 



BD+69 0993 


B 


10.2 


1.0 


G 








BD+08 4386 


A 


10.9 


0.7 


GO 








BD+08 4386 


B 


10.9 


0.4 


GO 








BD+03 4428 


A 


9.99 


0.77 


G5 








BD+03 4428 


B 


10.20 


0.80 


G5 








BD+76 0835 


A 


9.98 


0.81 


G 








BD+76 0835 


B 


10.04 


0.64 


GO 








BD+01 4575 


A 


10.7 


0.6 


GO 








BD+01 4575 


B 


10.7 


0.6 


GO 








BD+11 5033 


A 


10.26 


0.87 


KO 


8.99 


0.37 


0.47 


BD+11 5033 


B 


10.11 


0.93 


KO 


8.91 


0.29 


0.35 


HD 224984 




8.50 


0.55 


F8 








HD 224994 




8.72 


0.56 


GO 









References. — (CS99) Cuypers & Seggcwiss 1999; (D84) Duncan 1984; (M76) Mechler 1976; (MM78) 
Mechler & McGinnis 1978; (094) Olscn 1994 ; (S93) Schuster, Parrao & Contreras Martinez 1993 
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Table 3. Temperatures and Lithium Abundances of Program Stars. 



HD/BD 


T cff (J- K) 


T eff (fit) 


[m/H] 


EW (Li I) 


logN(Li)coG 


logN(Li) flt 




(K) 


(K) 




mA 






HD 4552 A 


6140 


6000 


-0.2 


<2 


<1.7 


<1.6 


BD+12 0090 


5837 


5750 


-0.2 


99+3 


2.2 


2.3 


HD 6872 A 


6543 


6250 


-0.2 


<3 


<2.0 


<1.8 


HD 6872 B 


5954 


6250 


-0.2 


36+3 


2.4 


2.6 


HD 8624 




5250 


-0.2 


31+4 




1.6 


HD 8610 




5500 


-0.2 


17+3 




1.6 


BD+60 269 




6250 


0.0 


49+4 




2.6 


BD+60 271 




6250 


0.0 


<3 




<1.8 


HD 31208 




5000 


-0.1 


<2 




<1.0 


BD+07 754 




5000 


-0.1 


<2 




<1.0 


HD 39274 


5836 


6000 


0.0 


61+3 


2.4 


2.6 


HD 39275 


5305 


5750 


-0.1 


79+3 


2.1 


2.5 


HD 54046 




6000 


-0.5 


38+3 




2.4 


HD 54100 




6250 


-0.5 


39+3 




2.6 


BD+28 1698 




6000 


-0.2 


65+2 




2.7 


BD+28 1697 




6000 


-0.2 


67+2 




2.7 


BD+15 2080 




5750 


-0.2 


<2 




<0.7 


BD+15 2079 




6000 


-0.2 


9+2 




1.7 


BD+34 2091 A 


5723 


5750 


-0.1 


88+4 


2.6 


2.6 


BD+34 2091 B 


5954 


5750 


-0.1 


63+3 


2.6 


2.5 


HD 92222 A 


5257 


5500 


-0.1 


24+3 


1.4 


1.7 


HD 92222 B 


5404 


5500 


0.0 


8+2 


1.4 


1.4 


BD+13 2311 A 


6403 


6500 


0.0 


37+3 


2.8 


2.8 


BD+13 2311 B 


6403 


6500 


0.0 


<2 


<1.9 


<2.0 


HD 98744 


6077 


6250 


-0.2 


<2 


<1.5 


<1.5 


HD 98745 


6140 


6250 


-0.2 


36+3 


2.4 


2.5 


HD 111484 A 




6000 


-0.2 


64+3 




2.6 


HD 111484 B 




6000 


-0.2 


92+3 




2.8 


HD 124257 A 


6139 


6000 


0.0 


65+3 


2.7 


2.6 


HD 124257 B 


5954 


5750 


0.0 


11+3 


1.9 


1.7 


HD 124913 A 


6472 


6500 


0.0 


<4 


<2.0 


<2.2 


HD 124913 B 


6763 


6500 


-0.1 


<4 


<2.3 


<2.1 


BD+02 2820 A 




6000 


0.0 


<3 




<1.5 


BD+02 2820 B 




5500 


0.0 


<3 




<1.2 


BD+35 2576 A 




6000 


0.0 


8+2 




1.8 


BD+35 2576 B 




6000 


0.0 


4+2 




1.7 


BD+13 2830 A 


5668 


5750 


-0.2 


<5 


<1.3 


<1.7 


BD+13 2830 B 


5354 


5750 


-0.2 


<3 


<1.0 


<1.4 


BD+15 2867 B 


5210 


5250 


-0.2 


<3 


<1.0 


<1.0 


BD+15 2867 C 


5070 


5250 


-0.2 


<2 


<0.8 


<0.9 


HD 155674 




5250 


-0.2 


<5 




<1.1 


BD+54 1862 




5250 


-0.2 


<3 




<1.1 


BD+65 1043 


4615 


5250 


0.0 


<2 


<0.4 


<1.2 


BD+65 1044 


6688 


6500 


0.0 


<3 


<1.6 


<1.5 


BD+74 0718 




5500 


-0.1 


<3 




<1.2 


BD+74 0719 




5500 


-0.1 


<3 




<1.2 


HD 167215 




6500 


-0.2 


43+3 




2.8 


HD 167216 




6500 


-0.2 


34+3 




2.7 
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Table 3 — Continued 



HD/BD T cH (J-K) T eff (fit) [m/H] EW (Li I) logN(Li) CO G logN(Li) flt 



BD+69 0993 A 




6500 


0.0 


30±3 




2.7 


BD+69 0993 B 




6500 


0.0 


34±3 




2.7 


BD+08 4386 A 




6000 


0.0 


<2 




<1.6 


BD+08 4386 B 




5500 


0.0 


<3 




<1.2 


BD+03 4428 A 




5500 


-0.2 


<3 




<1.2 


BD+03 4428 B 




5500 


-0.2 


<3 




<1.4 


BD+76 0835 A 




6000 


-0.2 


<2 




<1.5 


BD+76 0835 B 




5750 


-0.2 


<2 




<1.3 


BD+01 4575 A 




5750 


-0.2 


<3 




<1.2 


BD+01 4575 B 




5750 


-0.2 


<3 




<1.6 


BD+11 5033 A 


5162 


5250 


-0.1 


<3 


<1.0 


<1.2 


BD+11 5033 B 


5779 


5000 


-0.1 


<3 


<1.4 


<1.1 


HD 224984 




6000 


-0.1 


35±3 




2.3 


HD 224994 




6000 


-0.1 


33±3 




2.3 
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ABSTRACT 

We present high-resolution spectroscopic observations of the Li I resonance 
line in a sample of 62 stars that belong to 31 common-proper motion pairs with 
twin F or G-type components. Photospheric abundances of lithium were derived 
by spectral synthesis analysis. For seven of the pairs, we have measured large 
lithium abundance differences. Eleven other pairs have components with similar 
lithium abundances. We cannot determine if the remaining 13 pairs have lithium 
differences because we did not detect the Li I lines, and hence we can only provide 
upper limits to the abundances of both stars. Our results demonstrate that twin 
stars do not always share the same lithium abundances. Lithium depletion in 
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solar-type stars does not only depend on age, mass, and metallicity. This result 
is consistent with the spread in lithium abundances among solar-type stars in 
the solar-age open cluster M67. Our stars are brighter than the M67 members 
of similar spectral type, making them good targets for detailed follow-up studies 
that could shed light on the elusive mechanism responsible for the depletion of 
lithium during the main-sequence evolution of the Sun and solar-type stars. 

Subject headings: stars: abundances — stars: late-type - binaries: visual - stars: 
evolution - stars: fundamental parameters 

1. Introduction 

Lithium has a rich life. It is created in many environments (primordial nucleosynthe- 
sis, cosmic ray spallation in the interstellar medium, spallation and fusion reactions around 
compact relativistic objects, flares, and red giant thermal pulses), and it is destroyed in the 
interior of stars by collision with protons at a temperature of about 2.7 xlO 6 K. The Sun 
has a lithium abundance more than two orders of magnitude lower than meteoritic material 
(Miiller, Peytremann, & de la Reza 1975), which reflects the composition of the presolar neb- 
ula. The Sun probably did not deplete lithium during the pre main-sequence evolution, but 
rather during the slow main-sequence evolution, although the exact mechanism has not been 
identified yet (Martin 1997, 1998). Recent models of lithium depletion in solar-type stars 
suggest that two types of mixing may be at work during main-sequence evolution; namely 
overshooting and rotational mixing (Umezu & Saio 2000). The history of rotational mixing 
may be strongly influenced by the initial conditions of angular momentum distribution in 
the protoplanetary disk. Thus, it is conceivable that there could be a connection between 
the rate of lithium depletion by rotational mixing and the presence of companions (stellar 
or substellar) to the stars. 

Binary stars with twin components are interesting for understanding binary formation 
and evolution. Common proper motion (CPM) pairs of twins allow us a test of the validity of 
stellar evolution models because their separations are so large (>100 AU) that each member 
has probably evolved independently. If all the properties of stars are determined by age, 
mass, and chemical composition, twins should be identical. If the lithium abundance of 
a main-sequence star is determined by its age, mass, and chemical composition, pairs of 
twins should show the same lithium abundances. In the early work of Herbig (1965), it 
was already apparent that there was a scatter in lithium abundances among pairs of twins. 
Herbig's study included 53 UMa A and B, where lithium was detected in A but not in B; 
16 Cyg A and B, where lithium was detected in A but not in B; and 53 Aqr A and B, where 
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lithium was detected in both stars. 

The system 16 Cyg has received some attention recently. The two stars of the pair are 
nearly identical. The primary has V=5.96, T e g=5785 K and log g=4.28. The secondary has 
V=6.20, T e ff=5747 K and log g=4.35. Despite those similarities, the members of this pair 
have two fundamental differences. King et al. (1997) have shown that the primary has a 
lithium abundance higher than that of the Sun (logN(Li)=1.27±0.05 in the customary scale 
of log N(H)=12), while the secondary has a much lower lithium abundance (logN(Li)<0.60). 
The large difference in lithium abundance is surprising because both stars should be the same 
age and have very similar masses. Another significant difference between the two stars is that 
the secondary harbors a giant planet with a minimum mass of 1.5 Jupiters, an eccentricity 
of e=0.63 and a period of 800.8 days (Cochran et al. 1997). On the other hand, no giant 
planet has yet been detected around 16 Cyg A. However, such nondetection does not rule 
out the presence of giant planets with long periods or smaller planets with short periods. 

Gonzalez & Laws (2000) and Ryan (2000) have studied the distribution of lithium 
abundaces among planet-harboring stars, and have reached opposite conclusions. While 
Gonzalez & Laws claim that lithium is overdepleted in stars with giant planets, Ryan claims 
that it is normal. The disagreement comes from the choice of stars to be compared with the 
host stars of extrasolar planets. This is a difficult task because the ages of most field stars 
are not well known. 

Pairs of twins offer the possibility of comparing stars that are coeval, and have nearly 
the same mass and metallicity. Thus, we may expect that the systematic study of binaries 
may shed light on whether there is any connection between the presence of giant planets 
and lithium depletion. Gratton et al. (2001) have recently reported an abundance analysis 
for six wide binaries. They find that the primary star of the binary HD 219542 has higher 
iron and lithium abundances than the secondary star, which leads them to suggest that the 
primary has ingested a planet. 

In this paper we report lithium abundances for a sample of 31 visual binaries. The tar- 
gets have been chosen to be as similar as possible to 16 Cyg in spectral type and magnitude. 
The paper is organized as follows: In Section 2 we describe the observations. In Section 3 we 
present the abundance analysis. In Section 4 we discuss each pair of twins individually when 
we could detect the Li I line in at least one of the stars, and we compare them with lithium 
abundances in open clusters and theoretical models. Section 5 contains the discussion, and 
in Section 6 we present our conclusions. 
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2. Observations 

We selected our targets from the list of CPM pairs compiled by Halbwachs (1986). We 
filtered the sample by requiring that the difference in V magnitudes should be less than 0.4, 
and the spectral type should be between F8 and K0. The declination of the stars had to be 
north of —30 degrees. Our sample is listed in Table 1. 

Observations were obtained during several observing runs at Lick and Keck observa- 
tories between January 1998 and March 2001. The binaries BD+08 4386, BD+69 0993, 
BD+74 0718, and BD+76 0835 were observed at Keck. All the other binaries were observed 
at Lick. The observing log is given in Table 1. 

The Lick observations were carried out using the Hamilton echelle spectrograph on the 
Shane 3-meter telescope. This instrument provides 92 orders covering the spectral range 
from 490 nm to 890 nm at a resolution of 48,000. A standard reduction was performed 
(Valenti 1994), which included background subtraction, cosmic ray removal, and flat fielding. 
Wavelength calibration was made using ThAr lamp exposures. 

The Keck observations were obtained with the HIRES echelle spectrograph (Vogt et al. 
1994). The resolution was 67,000 and the coverage was from 600 to 850 nm. Data reduction 
was performed using IDL routines similar to those used for the Lick data. 

3. Analysis 

We have compiled the available photometry in the literature for our program stars and 
we provide it in Table 2. The V magnitudes, B — V colors, and spectral types were taken from 
SIMBAD, except for the particular cases listed in Section IV, and for the notes to Table 2. 
This photometric database is rather heterogeneous. Some of our stars are included in the 
2MASS second incremental release point source catalog (PSC). Their 2MASS J — K colors 
are included in Table 2. Temperatures were derived using the empirical scale of Alonso, 
Arribas & Martinez-Roger (1996) for solar metallicity. The temperatures obtained from the 
J — K colors are given in Table 3. Typical error bars in J — K colors are ±0.03, corresponding 
to ±180 K in T e g. These temperatures were derived to compare with our, presumably more 
reliable, temperatures obtained from synthetic spectra fitting to the observed spectra in the 
lithium region. This comparison is shown in Figure 1 for 26 stars that have available 2MASS 
JK photometry. The mean difference between the T e g derived from J — K and from fits 
is —65 K and the standard deviation is 270 K. We conclude that there is no significant 
systematic deviation of the T e g derived from J — K colors with respect to those derived from 
synthetic fits. 
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For each spectrum we determined the continuum level using the DECH20 software 
(Galasutdinov 1992). To detect the real continuum level, we smoothed the observed spectra 
using three loops of a 3-point smoothing boxcar. Then the continuum points were fitted 
by splines to use in the fine analysis procedure. The S/N ratios per pixel given in Table 1 
were obtained by measuring the standard deviation of the normalized continuum in the 
spectral range 669.3-669.5 nm. We carried out the computations using the LTE spectral 
synthesis program WITA6 (Pavlenko 2000), which is a modified version of the program 
used by Pavlenko et al. (1995). Our computations of synthetic spectra were carried out 
using plane-parallel model atmospheres in LTE, with no energy divergence. The model 
atmospheres were taken from Kurucz's grid available on the web. 

We used metallicities in the range [m/H]=0.0 to [m/H]=-0.2 for fitting the observed 
spectra 1 except for HD 54046 and HD 54100, which are known metal poor stars (Norris 
1986). Chemical equilibrium was computed for the mix of RilOO molecular species, with 
data taken from Tsuji (1973). ATLAS9 (Kurucz 1993) grid of continuum opacity sources 
were used in our computations. We used VALD (Piskunov et al. 1995) for the atomic line 
list and damping constants. For some lines the damping data are missing; in that case we 
used Unsold (1955) approach to compute them. We had to increase the log gf of the Fel 
line at 670.510 nm by a factor of 2 (from 0.0319 given by VALD to 0.0640) to find good 
fits to the observed spectra of the Sun as a star (Kurucz et al. 1984). Computations were 
carried out for a fixed microturbulent velocity of V t =2 km s _1 and a Voigt profile for every 
absorption line. 

Synthetic spectra were computed with a step of 0.002 nm. Instrumental profile and 
macroturbulence broadening effects were modeled by convolution with a Gaussian of half- 
width in the range 0.010 nm to 0.018 nm. For most stars we could desribe the profile 
with instrumental broadening only, without any rotation. That provides upper limit of 
vsim=3 km s -1 . No enhanced van-der Waals broadening was used. For two stars rotational 
broadening was taken into account (HD 6872 A and BD+65 1044; vsim=14 km s _1 ). For 
each star we found the best-fitting synthetic spectrum using a least-squares minimization 
algorithm. Figures 2, 3 and 4 show typical examples of synthetic fits to the echelle spectra. 

For each LTE lithium abundance, we determined NLTE abundance corrections usind 
the grid of curves of growth of Pavlenko & Magazzu (1996). In Table 3 we summarize our 
results. We give the temperatures, metallicities and lithium abundances (including NLTE 
corrections) that we used in the best fits to the observed spectra. Equivalent widths of the Li I 
line measured by direct integration of the line profile in the observed spectrum are given in 



1 We refer to the solar abundances given by Grevesse (1984) using the convention [m/H]=0 
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Table 3 as well. They can be used to derive lithium abundances from curve of growth (COG) 
computations. For the stars with 2MASS photometry we have obtained COG Li abundances 
using the temperatures derived from the (J — K) — T e g calibration of Alonso et al. (1996) and 
the WITA code. Those abundances are listed in Table 3 for comparison with the abundances 
obtained from the synthetic spectra. Figure 5 illustrates this comparison for the same 26 
stars that were used in Figure 1. The mean difference between the lithium abundances 
derived using these two methods is logN \Li) COG — log]S {Li)^ = —0.08 dex, and the standard 
deviation is 0.17 dex. We conclude that there is no significant systematic discrepancy between 
the two methods of obtaining lithium abundances. We adopt the uncertainty of 1 ct=0.17 dex 
in our absolute lithium abundances. In our computations we have assumed that our stars 
are dwarfs. A variation of log g ±0.5 dex provides a maximum AlogN(Li)<0.1 dex. 



4. Comments on Individual Stars 

The main goal of this paper is to determine whether stars with the same age and initial 
chemical composition, and nearly equal masses, have depleted lithium at the same rate. We 
limit our discussion to pair of twins for which we could detect the Li I line in at least one of 
the stars. When we could not detect the Li I line in any of the stars, we could not tell if the 
lithium abundances are different or similar, and therefore those pairs are not discussed any 
further. After a literature search using SIMBAD, we found that some of the pairs deserve 
individual comments: 

• HD 6872 A and B: Oblak & Chareton (1980) and Duncan (1984) included HD 6872 A 
and B among a sample of stars that are sufficiently evolved for the ages to be estimated 
from T e g and absolute magnitudes using the displacement from the zero-age main- 
sequence (ZAMS). Both groups obtained uvby/3 photometry and derived T e g of 6380 K 
(Oblak & Chareton) and 6300 K (Duncan) for the primary. For the secondary, they 
gave 6251 K (Oblak & Chareton) and 6084 K (Duncan). The discrepancy in T e g 
is mainly due to their use of different equations of b — y index versus temperature. 
Soderblom, Duncan & Johnson (1991) converted the b — y color of Duncan to B — V 
color. In Table 1, we adopt the Soderblom et al. color, which is different from the 
color provided by SIMBAD. Our best fit temperatures are 6250 K for both stars. We 
find good synthetic fits to our spectra with the same T e g for both stars (Figure 3). 
The fainter star has a higher lithium abundance than the hotter star by at least 0.8 
dex. Adopting the T e ff of either Oblak & Chareton or Duncan does not change the 
conclusion that the lithium abundances of these two stars are significantly different. 
The primary star has broader lines than the secondary. Our best fit to the spectrum 



- 7- 



of HD 6872 A has rotational broadening of vsim=14 km s 1 . 

• HD 39274 and HD 39275: The lithium abundance difference between these two stars 
is 0.1 dex, which is less than the 1 a uncertainty in our analysis. HD 39274 is a binary 
with separation 0.225 arcsec (McAlister et al. 1989). 

• HD 54046 and HD 54100: Metal poor stars. Norris (1986) reported [Fe/H]=-0.6 and 
Eggen (1998) gave [Fe/H]=-0.45 based on photometric colors. 

• HD 8624 and HD 8610: Tokovinin (1999) found that HD 8624 is a double-lined spec- 
troscopic binary with a period of 14.9 days. HD 8610 has a K-type companion, and 
therefore the system is a quadruple. Tokovinin gave V=8.83 for both components of 
the HD 8624 close system, and inferred a spectral type of G3. 

• HD 98744 and HD 98745: HD 98745 has a lithium abundance 1 dex or more higher 
than HD 98744. Duncan (1984) found T eff =6050 K and 6039 K for HD 98744 and HD 
98745, respectively. Our best fit T e g is 6250 K for both stars. Adopting the T e g of 
Duncan does not change the conclusion that the lithium abundances of these two stars 
are significantly different. 

• HD 167215 and HD 167216: In Table 2 we give the B - V colors given by SIMBAD, 
because they are more consistent with the temperatures that we find in our analysis, 
than the B — V colors of Soderblom et al. (1991). These latter authors derive T e ff of 
5984 K and 6040 K for HD 167215 and HD 167216, respectively. We find that we need 
hotter temperatures to fit the spectra in the lithium region. We also need subsolar 
metallicities. The lithium abundance difference is 0.1 dex. The Li I resonance line 
of HD 167215 is clearly stronger than that of HD 167216. According to Soderblom 
et al., HD 167215 is more evolved above the ZAMS than HD 167216. The gravity of 
HD 167215 may be somewhat lower than that of HD 167216, and this could explain 
the difference in Li I line strengths. 

• HD 224984 and HD 224994: These two stars have nearly identical spectra. HD 224994 
itself is a close binary (separation=0.24 arcsec in 1988; McAlister, Hartkopf, & Franz 
1990), so the system is at least triple. The two stars have very similar spectra, including 
the Li I lines. We obtain lithium abundances of logN(Li)=2.3, which is about 0.8 dex 
below the cosmic value. Therefore, this is an example of a pair of twins that have 
depleted the same amount of lithium. The presence of a low-mass star companion at 
about 20 AU of HD 224994 has not changed the lithium depletion rate of this star with 
respect to HD 224984. 



- 8 - 



5. Discussion 

Lithium abundance is very sensitive to the mass (and thus T e ff) of stars. In the temper- 
ature range of our program stars (6500-5250 K), lithium abundances generally decrease with 
decreasing T e g, although at the hot end of this temperature range, the F-type lithium dip 
starts to kick in (Boesgaard & Tripicco 1987). If the lithium abundances of the components 
of pairs of twins were dominated by the mass dependance of lithium depletion, we would 
expect a correlation between the difference in T e g and the difference in logN(Li), in the sense 
that the cooler components should show systematically lower lithium abundances than the 
hotter ones. Figure 6 shows that such correlation is not present in the data. We conclude 
that the observed pattern of lithium abundances in pairs of twins cannot be explained solely 
with the dependance of lithium depletion on stellar mass. 

We have found five pairs with nearly identical T e ff and nearly identical logN(Li), and 
four pairs with nearly identical T e g and very different logN(Li). Hereafter we call these 
latter four pairs "16 Cyg analogs", namely: HD 6872 A/B, BD+60 0269/BD+60 0271, 
BD+13 2311 A/B, and HD 98744/HD 98745. The common characteristics of "16 Cyg analogs' 
is that they are pairs of twins with lithium abundance difference between the two stars of the 
pair that exceeds 0.5 dex (which corresponds to 3 a significance in our synthetic fit analysis) 
and difference in T e g <200 K. 

Figure 7 illustrates the comparison between the pairs of twins and the Hyades cluster in 
a T e ff versus logN(Li) diagram. The Hyades single star members define a tight relationship 
between T e g and logN(Li). In general the twins do not follow the Hyades relation. A possible 
explanation for this discrepancy could be that there are many tidally locked binaries (TLB) 
hidden among the twins. We do not know of any example of a TLB in our sample. 

HD 8610 has a lithium abundance similar to Hyades members with the same T e ff. The 
lithium-T e ff locus in the Hyades declines steeply from 5500 K to 5250 K, so one would expect 
that HD 8624 should have a lithium abundance much lower than that of HD 8610, which 
is contrary to what we have found. HD 8624 is a double-lined spectroscopic binary with a 
period of 14.91 days (Tokovinin 1999) and an eccentricity e=0132. Even though it is not 
strictly a TLB because the orbit is not completely circularized, it may have preserved lithium 
in a manner similar to the TLBs studied by Barrado y Navascues et al. (1997), which show 
a trend of higher lithium abundances when compared to single stars. However, TLBs are 
relatively rare among the general population of solar-type stars. A program to monitor the 
radial velocity of 10 of our pairs has been started using HIRES with the iodine cell. There 
are already enough data for 16 of the stars to find TLBs, but none has been seen (Marcy 
2002, private communication). We conclude that it is very unlikely that TLBs can account 
for most of the discrepancy between the pairs of twins and the Hyades. 
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The average age of the pairs of twins is probably larger than the Hyades cluster age 
(^600 Myr). It is informative to compare the twins with the well-studied solar-age cluster 
M67. In Figure 7 we compare the pairs of twins with the members of M67. The M67 stars 
show a large spread of lithium abundances at a given T e g, which is similar to the behavior 
exhibited by our binaries. The evolution of lithium from the Hyades to M67 implies a 
mechanism of lithium depletion during the main-sequence lifetime of solar-type stars that 
is not very mass dependent. We have not found in the literature any models that can 
explain the evolutionary pattern of lithium abundances implied by the comparison between 
the Hyades and the M67 solar-type stars. We conclude that there is currently not a good 
theoretical explanation for the lithium abundances observed in our pairs of twins and the 
open clusters. In future papers we plan to tackle this problem observationally via searches 
for connections between the lithium abundances in the pairs of twins and the presence of 
companions (stellar or substellar), chromospheric activity, rotation, and the abundances of 
other chemical elements. At the moment it would be premature to suppose that the presence 
of planets had anything to do with the spread in lithium abundances. 

6. Final Remarks 

We report high-resolution spectroscopy observations of 62 solar- type stars in 31 common 
proper motion pairs of twins. These binaries are sufficiently separated that the evolution 
of each star is not expected to be influenced by its stellar companion. Lithium abundances 
have been derived from theoretical fits to the observed spectra using model atmospheres. 
A comparison with COG abundances derived from Li I equivalent widths and temperatures 
obtained from near-infrared colors shows that the absolute lithium abundances are reliable 
to within 0.17 dex (1 a). 

We have found seven pairs of twins with large lithium abundance differences. The pat- 
tern of lithium abundances in our sample of twins seems to be midway between those of the 
Hyades (tight lithium-T e fj relation) and M67 (large lithium spread for any T e g) open clus- 
ters. We have not found any evolutionary models in the literature that provide a satisfactory 
explanation for the observed spread in lithium abundances in solar-type stars. We plan to 
follow up with additional studies of our sample of twins to obtain an overall understanding 
of all the effects that may come into play in the evolution of solar-type stars. 

This paper is based on data collected with the 3-m Shane telescope at Lick Observatory, 
California, and the Keck I telescope at the W. M. Keck Observatory, Mauna Kea, Hawaii. 
The Keck observatory is operated as a scientific partnership among the California Institute 
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of Technology, the University of California, and the National Aeronautics and Space Admin- 
istration. The Observatory was made possible by the generous financial support of the W. 
M. Keck Foundation. This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. The authors wish to extend special thanks to those of Hawaiian 
ancestry on whose sacred mountain of Mauna Kea we are privileged to be guests. Without 
their generous hospitality, the Keck telescope observations presented therein would not have 
been possible. We thank Louise Good for correcting the English of the manuscript, and 
Geoff Marcy for communicating results of the radial velocity program for a subset of our 
sample. We are also indebted to Ann Boesgaard and David Barrado y Navascues for sending 
data files for lithium abundances of M67 stars, and Nina Kharchenko for astrometric data. 
The project was partly supported by a SRG grant of AAS to YP. 
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Fig. 1. — Comparison of effective temperatures obtained from J — K colors and from the 
best synthetic fits to the observed spectra. 
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Fig. 2. — Solid lines: Lick Hamilton echelle spectra of the CPM pair HD 4552 A and B. 
Dashed lines: Synthetic spectra for T e g=6000 K, [m/H]=— 0.2, and logN(Li)=1.6 (above), 
and for T eff =5750 K, [m/H]=-0.2, and logN(Li)=2.3 (below). 
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Fig. 3. — Solid lines: Lick Hamilton echelle spectra of binary HD 6872 A and B. 
Other name for HD 6872 B is BD+12 0090. Dashed lines: Synthetic spectra for 
T eff =6250 K, [m/H]=-0.2, and logN(Li)=1.8 (above), and for T eff =6250 K, [m/H]=-0.2, 
and logN(Li)=2.6 (below). 
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671.5 



Fig. 4. — Solid lines: Lick Hamilton echelle spectra of the CPM pair HD 9518 A 
(BD+60 0269) and B (BD+60 0271). Dashed lines: Synthetic spectra for T eff =6250 K, 
[m/H]=0.0, and logN(Li)=2.55 (above), and for T eff =6250 K, [m/H]=0.0, and logN(Li)=1.8 
(below). 
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Fig. 5. — Comparison of lithium abundances obtained with J — K colors and COG calcula- 
tions and with the best synthetic fits to the observed spectra. 
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Fig. 6. — Differences in the T e g adopted for pairs of twins versus the differences in the lithium 
abundances derived by us (solid symbols) and by previous investigations (open symbols). 
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Fig. 7. — Lithium abundances versus T e g for pairs of twins. Our lithium measurements are 
given as solid dots, with the components of each binary linked by a solid line. The effective 
temperatures of some pairs are shifted by 30 K to avoid confusion. The components of 16 Cyg 
are denoted with an open dot and a solid line. The locus of lithium abundances in single 
stars members of the Hyades cluster is shown with a dotted line. The Hyades data have been 
taken from Thorburn et al. (1993). Lithium abundances for M67 members are denoted with 
open hexagons for detected Li I lines, and asterisks for upper limits. Data for M67 stars come 
from Hobbs & Pilachowski (1986), Spite et al. (1987), Garcia Lopez, Rebolo, & Beckman 
(1988), Deliyannis et al. (1994), Pasquini, Randich, & Pallavicini (1997), Jones, Fischer, & 
Soderblom (1999), Barrett et al. (2001), and Barrado y Navascues, & Balachandran (2002). 
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Table 1. Observing Log. 



HD/BD 


Telescope 


Date 


Texp (s) 


S/N 


T-JTl A K KO 

riJJ 40 Oz 


o-m ijicK 


1998 


Nov 23 


1 Knn 
10UU 


i A n 
14U 


dlj-\- iz uuyu 


3-m Lick 


1998 


Nov 23 


loUU 


1UO 


Ml J Oo /z A 


o-m lick 


1998 


Sept 4 


1 nnn 
1UUU 


i A n 
14U 


XlU Do i Z Jj 


o-m J^ICK 


1998 


Sept 4 


1 1UU 


i Qn 
loU 


rliJ oOZ4 


o-m ijicK 


1998 


Sept 3 


DUU 


int; 
1UO 


TX"Pi Qcm 
riJJ oOlU 


o-m IjicK 


1998 


Sept 3 


ouu 


i Qn 
loU 


JdJJH-dU zby 


o-m ijicK 


1998 


Sept 3 


oUU 


1 1 n 
11U 




3-m Lick 


1998 


Sept 3 


1 /inn 

14UU 


i nn 
1UU 


riJJ olzuo 


3-m Lick 


2001 


March 8 


1 onn 
loUU 


i f^n 
10U 


LSJJ + UV (04 


O-m ijlCK 


2001 


March 8 


"i onn 
loUU 


i nn 

iyu 


TX"Pl Q007/I 


o-m IjicK 


2001 


March 8 


o f^nn 
zOUU 


i nn 

iyu 


riJJ oyz I 


o-m IjicK 


2001 


March 8 


oonn 
zzUU 


I/O 


rl 1 ) 04U40 


o-m ijicK 


1998 


Jan 22 


loUU 


oon 

ZoU 


TTTi pr A -1 fin 

riJJ 041UU 


o-m IjicK 


1998 


Jan 22 


1 onn 
loUU 


q A n 
o4U 


\)\ )-\- Zo loyo 


3-m Lick 


1998 


Jan 21 


1 onn 
IzUU 


01 n 
Z1U 




o-m lick 


1998 


Jan 21 


1 onn 
IzUU 


onn 
zuu 


pn i i c onon 
r>JJ-r-10 zuou 


o-m lick 


1998 


Jan 22 


oa nn 
z4UU 


01 n 
zlU 


r>JJ-|-10 zU / y 


o-m ijicK 


1998 


Jan 22 


i onn 

loUU 


zOU 


rSJJ+o4 zUyl A 


o-m IjicK 


2001 


March 9 


onnn 
oUUU 


nn 

yu 


\)\ )-\- zuyi i) 


o-m ijicK 


2001 


March 9 


Q^nn 


n k 

yo 


LI 1 ~» nOOOO A 
rl 1 ) \}ZZZZ A 


o-m ijicK 


2001 


March 8 


Q^nn 

oDUU 


i a n 

14U 


TJT\ 00000 T3 


o-m IjicK 


2001 


March 8 


oOUU 


1 QK 

loO 


I) 1 ~» 1 1 Q OQ1 1 A 

dLJ-\- lo Zoll A 


o-m ijicK 


2001 


March 9 


k -i on 

OloU 


OO K 

zzo 


PFt 1 1 Q OQ1 1 "D 


o-m IjicK 


2001 


Jan 22 


A ACiCi 
44UU 


oon 
zzU 


riJJ y<5 i 44 


o-m IjicK 


2001 


March 8 


A nnn 
4UUU 


1(0 


II 1 ~» OQ7/1C; 

n 1 ) yo / 40 


o-m ijicK 


2001 


March 8 


a Knn 
40UU 


1(0 


TJ"Pl 1 1 1 A Q A A 
riJJ 1114o4 A 


o-m IjicK 


1998 


Jan 21 


1 onn 
loUU 


onn 
zUU 


II 1 ~» 1 1 1 A OA L) 

rl 1 ) 111404 tS 


o-m ijicK 


1998 


Jan 21 


i onn 

loUU 


ocn 

ZoU 


II 1 ~» 10/10K7 A 


o-m ijicK 


1999 


May 26 


oa nn 

z4UU 


i on 

loU 


riJJ Iz4z0( ri 


o-m IjicK 


1999 


May 26 


1 onn 
loUU 


i c;n 
10U 


LI 1 ~» i O/l Ol Q A 
1IU lZ4yio A 


o-m ijicK 


1999 


May 27 


i onn 

loUU 


1 1 n 
11U 


riJJ lz4yio ri 


o-m IjicK 


1999 


May 27 


1 onn 
loUU 


1U0 


SdLJ-\-\jZ ZoZU a 


o-m ijicK 


1999 


May 26 


onnn 
zuuu 


Q K 

yo 


f» n _i_ n 9 9»9n pi 

DIJyUZ ZoZU -D 


o-m IjIck 


1999 


May 26 




120 


BD+35 2576 A 


3-m Lick 


1999 


May 27 


2200 


135 


BD+35 2576 B 


3-m Lick 


1999 


May 27 


2000 


100 


BD+13 2830 A 


3-m Lick 


2001 


March 8 


2000 


65 


BD+13 2830 B 


3-m Lick 


2001 


March 8 


2400 


60 


BD+15 2867 B 


3-m Lick 


1999 


May 28 


2200 


80 


BD+15 2867 C 


3-m Lick 


1999 


May 28 


2400 


85 


HD 155674 


3-m Lick 


1999 


May 28 


2000 


85 


BD+54 1862 


3-m Lick 


1999 


May 28 


2000 


110 


BD+65 1043 


3-m Lick 


1998 


Sept 4 


1800 


105 


BD+65 1044 


3-m Lick 


1998 


Sept 4 


2000 


100 


BD+74 0718 


10-m Keck I 


1998 


Sept 20 


200 


95 


BD+74 0719 


10-m Keck I 


1998 


Sept 20 


200 


150 


HD 167215 


3-m Lick 


1998 


Sept 3 


800 


120 


HD 167216 


3-m Lick 


1998 


Sept 3 


800 


110 


BD+69 0993 A 


10-m Keck I 


1998 


Sept 20 


120 


155 
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Table 1 — Continued 



HD/BD 




Telescope 


Date 




Texp (s) 


S/N 


BD+69 0993 


B 


10-m Keck I 


1998 


Sept 


20 


120 


115 


BD+08 4386 


A 


10-m Keck I 


1998 


Sept 


20 


180 


105 


BD+08 4386 


B 


10-m Keck I 


1998 


Sept 


20 


180 


110 


BD+03 4428 


A 


3-m Lick 


1998 


Sept 


4 


1500 


120 


BD+03 4428 


B 


3-m Lick 


1998 


Sept 


4 


1800 


110 


BD+76 0835 


A 


10-m Keck I 


1998 


Sept 


20 


180 


160 


BD+76 0835 


B 


10-m Keck I 


1998 


Sept 


20 


180 


150 


BD+01 4575 


A 


3-m Lick 


1998 


Sept 


4 


2500 


80 


BD+01 4575 


B 


3-m Lick 


1998 


Sept 


4 


2500 


80 


BD+11 5033 


A 


3-m Lick 


1998 


Sept 


3 


1200 


100 


BD+11 5033 


B 


3-m Lick 


1998 


Sept 


3 


1500 


105 


HD 224984 




3-m Lick 


1998 


Sept 


3 


800 


110 


HD 224994 




3-m Lick 


1998 


Sept 


3 


1000 


90 
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Table 2. Photometric Data for Program Stars. 



HD/BD 


V 


B - V 


Sp.T. 


J 


J - H 


J - K 


b-y 


ml 


cl 


Sources 


TTTl A EEO 

riJJ 400z 


O 07 

o.o / 


n eq 
U.Oo 


T?Q 

r o 


1 Q/1 
i .o4 


n oo 
U.zz 


n on 

u.zy 










■ > i ^ — t — i — uuyu 


o on 
y .zU 


n E£ 

U.OO 


*._tU 


O.Ul 


n Qn 
U.oU 


n Q/i 
U.o4 










TTTl C07O A 
riJJ Do f Z A 


Q nfi 
o.Ub 


n a 1 
U.41 


T?Q 

r o 


1 1 fi 

/ .lb 


n i fi 
U. lb 


n oq 
U.zo 


n om 

u.zy i 


n i fio 

u. ioy 


n a 01 
U.4zl 


A/T7fi 
1V1 ( 


TTTl fifi70 "R 
XlU Do (Z D 


fi EQ 
O.JO 


n Q7 

U.o f 


T?Q 

r o 


I .ou 


n oa 

U.Z4 


n qo 

U.oZ 


n qqk 

U.ooO 


n i ee 
u . 1 


n Q7K 

U.o I 


1V1 ( O 


un QfiO/i 


q m 

o.Ul 


n ^7 

U.O f 


VjU 








n /iO/i 

U.4Z4 


n OQ/i 

U.zo4 


n QQn 
U.ooU 




tt Ti ofii n 
liU oOlU 


q m 
o.Ul 


n fin 
U.OU 


T?Q 

r o 








n /inn 
U.4UU 


n i oi 

u. iy i 


n /i nE 
U.4U0 


uy4 


I) 1 ~» 1 fin Ofin 

riD-rbU zby 


O.Ob 


n /i q 
U.4o 


T?Q 

r o 
















id n i en 071 


o nc 
y .Uo 


n Ad 
U.4y 


T7 Q 
1 O 
















XT Tl QIOftQ 

liU olzUo 


q on 
o.zU 


n qq 
U.oo 


Ti'n 








n /i ofi 

u.4yo 


n QOQ 
U.oyo 


n q i o 
U.olz 


GOQ 


1) 1 ~» 1 ri'7 7E/1 

LSJJ + Ui /04 


Q Q Q 

o.oo 


n on 
U.oy 


rVU 








n E1 e 
U.Olo 


n /i O/i 
U.4z4 


n Qn7 
U.oU / 


COQ 


TTTl QQ07/I 


Q Ofi 

o.yo 


n ei 
U.01 


VjU 


7 QO 
i .o» 


n q 1 
U.ol 


n q a 
U.o4 










TTTl Qn07K 

liU oyz ( 


q on 

o.yy 


r\ qa 
U.o4 


VjU 


7 fiQ 
i .Oo 


n Qfi 
U.ob 


n A A 
U.44 










un E/10/1£ 


7 70 


U.D4 


CC\ 
\J\J 








n QQO 

U.ooy 


n i eq 
U. lOo 


n oo£ 
u.zyo 


TlQ/1 
Uo4 


tt ri E/i 1 nn 
riJJ 041UU 


1 70 


n kq 
U.oo 


T?Q 

r o 








n Qfio 
U.oOz 


n i eo 

u. ioy 


n q i fi 
U.olO 


JJo4 


r>D-|-zo loyo 


Q 07 
O.Z / 


n co 
U.oo 


VjrU 
















JJU-t-Zo luy i 


Q 0£ 
O.ZO 


U.OU 


*._tU 
















pn i i e onon 
rSJJH-lo zUoU 


o m 

y.ui 


n rq 
U.Oo 


/^i c- 
\jD 
















r>JJ-|-10 zU/y 


Q 77 


u.oy 


VjrU 
















rSJJ+o4 zUyi A 


O 71 

y . ( i 


n c;i 
U.01 


VjU 


q sn 
o.OU 


n on 

u.zy 


n Qfi 
U.oO 










I) 1 ~» 1 Q A onoi U 

i) i j -|- . > 4 zuyi i) 


O 70 

y. ( y 


U.4y 


VjU 


O.oO 


n oq 
U.zo 


n qo 
U.oz 










LI 1 ~» OOOO 1 ) A 

nu yzzzz a 


O AC\ 

y .4u 


n 7n 

U. i U 


^n 


Q 1 7 
o.l i 


n qq 
U.oo 


n /i e 
U.40 


n /i q q 

U.4oo 


n oqi 
U.zol 


n QOQ 

U.ozo 


TlQ/1 
Uo4 




o Qn 

y .ou 


n 7o 

u. i y 


VjU 


Q 1 /I 

o.l4 


n Qfi 
U.oO 


n /i o 
U.4z 


n /i Efi 
U.400 


n ofin 
U.zOU 


n Q7fi 
U.o i 


JJo4 


U>\ \ 1 1 Q OQ 1 1 A 

dLJ-\- lo Zoll A 


O QO 


U.4D 


T?Q 

r o 


Q 7/1 
o. ( 4 


n oa 

U.Z4 


n oe 
U.zO 










pn i iq oq i 1 ~n 
sdL)-\-io zoll ri 


Q Q 
O.O 


n 

U.o: 


T?Q 

r o 


Q Q K 

o.oO 


n oq 
U.zo 


n oe 
U.zO 










TJTl OQ7/I /I 
riJJ y<5 f 44 


y.io 


Ci Ad 

U.40 


VjU 


Q OO 

o.zz 


n ofi 
U.zO 


n Qn 
U.oU 


n QQfi 
U.ooO 


n i /i o 
U. 14z 


n /i i o 

u.4iy 


TVllVl i O 


riJJ ye ( 40 


n on 

y.oy 


n /i 7 
U.4 / 


VjrU 


Q /I O 

o.4z 


n oe 
U.zO 


n on 

u.zy 


n QQ E 

U.ooO 


n i eo 
U. loz 


O Q7Q 

U.o / o 


IVIJVI / O 


TJTl 1 1 1 A Q A A 
1 1 U 1114o4 A 


o 77 
O.I 1 




VjU 














Ooyy 


T-JTl 1 1 1 A Q/1 U 
rllJ 111404 1) 


Q QQ 
o.oo 




^n 

VjrU 
















TJTl 1 O/I OE7 A 
rl 1 J 1Z4Z0/ A 


O 1 

y . i 


n q. 
U.o. 


VjrO 


o in 
o.lU 


n oo 
U.zz 


n oo 
U.zy 










TJTl 10/10E7 "D 

riJJ Iz4z0i ri 


o fi 

y .0 


i 1 . 
l.l: 


\jD 


Q fiQ 

o.bo 


n oe 
U.zO 


n qq 
U.oo 










TJTl 1 O/I 01 Q A 
rl 1 ) lZ4yio A 


o n 

y .u 


l.U 


VjrO 


Q EO 

o.oy 


n 01 
U.Zl 


n O/i 

U.Z4 










T-JTl 10/101Q Ti 

\\\ ) lz4yio ri 


o o 

y .z 


1 1 
1.1 


VjO 


q no 

o.yz 


n i 7 
U. 1 / 


n on 
U.zU 










pn i no oQon A 
rSJJ+Uz ZozU A 


inn 
1U.U 


n n 

u.y 


VjO 
















L) i a i no OQon T3 
JJJJ+Uz zozU ri 


inn 
1U.U 


n n 

u.y 


*j0 
















BD+35 2576 A 


10.20 


0.67 


K0 
















BD+35 2576 B 


10.77 


0.61 


K0 
















BD+13 2830 A 


9.8 


1.1 


K0 


9.28 


0.34 


0.37 










BD+13 2830 B 


9.9 


1.1 


K0 


9.37 


0.33 


0.43 










BD+15 2867 B 


9.7 


1.1 


K0 


9.15 


0.39 


0.46 










BD+15 2867 C 


9.7 


1.1 


K0 


9.16 


0.40 


0.49 










HD 155674 


8.80 


1.20 


K0 
















BD+54 1862 


9.29 


1.24 


K8 
















BD+65 1043 


8.95 


0.99 


G5 


7.19 


0.46 


0.60 










BD+65 1044 


9.67 


0.43 


F8 


8.84 


0.18 


0.21 










BD+74 0718 


11.1: 


0.87 


GO 
















BD+74 0719 


10.09 


0.77 


G5 
















IID 167215 


8.10 


0.48 


F8 








0.357 


0.148 


0.373 


D84 


HD 167216 


8.23 


0.48 


F8 








0.346 


0.160 


0.357 


D84 


BD+69 0993 A 


10.1 


1.0 


G 
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Table 2 — Continued 



HD/BD V B - V Sp.T. J J - H J - K b-y ml cl Sources 



BD+69 0993 B 


10.2 


1.0 


G 








BD+08 4386 A 


10.9 


0.7 


GO 








BD+08 4386 B 


10.9 


0.4 


GO 








BD+03 4428 A 


9.99 


0.77 


G5 








BD+03 4428 B 


10.20 


0.80 


G5 








BD+76 0835 A 


9.98 


0.81 


G 








BD+76 0835 B 


10.04 


0.64 


GO 








BD+01 4575 A 


10.7 


0.6 


GO 








BD+01 4575 B 


10.7 


0.6 


GO 








BD+11 5033 A 


10.26 


0.87 


KO 


8.99 


0.37 


0.47 


BD+11 5033 B 


10.11 


0.93 


KO 


8.91 


0.29 


0.35 


HD 224984 


8.50 


0.55 


F8 








HD 224994 


8.72 


0.56 


GO 









References. — (CS99) Cuypers & Seggewiss 1999; (D84) Duncan 1984; (M76) Mechler 1976; (MM78) 
Mechler & McGinnis 1978; (094) Olsen 1994 ; (S93) Schuster, Parrao & Contreras Martinez 1993 



Table 3. Temperatures and Lithium Abundances of Program Stars. 



HD/BD 


T eff (J - K) 


T eff (fit) 


[m/H] 


EW (Li I) 


logN(Li)coo 


logN(Li) fit 




(K) 


(K) 




ml 






HD 4552 A 


6140 


6000 


-0.2 


<2 


<1.7 


<1.6 


BD+12 0090 


5837 


5750 


-0.2 


99+3 


2.2 


2.3 


HD 6872 A 


6543 


6250 


-0.2 


<3 


<2.0 


<1.8 


HD 6872 B 


5954 


6250 


-0.2 


36+3 


2.4 


2.6 


HD 8624 




5250 


-0.2 


31+4 




1.6 


HD 8610 




5500 


-0.2 


17+3 




1.6 


BD+60 269 




6250 


0.0 


49+4 




2.6 


BD+60 271 




6250 


0.0 


<3 




<1.8 


HD 31208 




5000 


-0.1 


<2 




<1.0 


BD+07 754 




5000 


-0.1 


<2 




<1.0 


HD 39274 


5836 


6000 


0.0 


61+3 


2.4 


2.6 


HD 39275 


5305 


5750 


-0.1 


79+3 


2.1 


2.5 


HD 54046 




6000 


-0.5 


38+3 




2.4 


HD 54100 




6250 


-0.5 


39+3 




2.6 


BD+28 1698 




6000 


-0.2 


65+2 




2.7 


BD+28 1697 




6000 


-0.2 


67+2 




2.7 


BD+15 2080 




5750 


-0.2 


<2 




<0.7 


BD+15 2079 




6000 


-0.2 


9+2 




1.7 


BD+34 2091 A 


5723 


5750 


-0.1 


88+4 


2.6 


2.6 


BD+34 2091 B 


5954 


5750 


-0.1 


63+3 


2.6 


2.5 


HD 92222 A 


5257 


5500 


-0.1 


24+3 


1.4 


1.7 


HD 92222 B 


5404 


5500 


0.0 


8+2 


1.4 


1.4 


BD+13 2311 A 


6403 


6500 


0.0 


37+3 


2.8 


2.8 


BD+13 2311 B 


6403 


6500 


0.0 


<2 


<1.9 


<2.0 


HD 98744 


6077 


6250 


-0.2 


<2 


<1.5 


<1.5 


HD 98745 


6140 


6250 


-0.2 


36+3 


2.4 


2.5 


HD 111484 A 




6000 


-0.2 


64+3 




2.6 


HD 111484 B 




6000 


-0.2 


92+3 




2.8 


HD 124257 A 


6139 


6000 


0.0 


65+3 


2.7 


2.6 


HD 124257 B 


5954 


5750 


0.0 


11+3 


1.9 


1.7 


HD 124913 A 


6472 


6500 


0.0 


<4 


<2.0 


<2.2 


HD 124913 B 


6763 


6500 


-0.1 


<4 


<2.3 


<2.1 


BD+02 2820 A 




6000 


0.0 


<3 




<1.5 


BD+02 2820 B 




5500 


0.0 


<3 




<1.2 


BD+35 2576 A 




6000 


0.0 


8+2 




1.8 


BD+35 2576 B 




6000 


0.0 


4+2 




1.7 


BD+13 2830 A 


5668 


5750 


-0.2 


<5 


<1.3 


<1.7 


BD+13 2830 B 


5354 


5750 


-0.2 


<3 


<1.0 


<1.4 


BD+15 2867 B 


5210 


5250 


-0.2 


<3 


<1.0 


<1.0 


BD+15 2867 C 


5070 


5250 


-0.2 


<2 


<0.8 


<0.9 


HD 155674 




5250 


-0.2 


<5 




<1.1 


BD+54 1862 




5250 


-0.2 


<3 




<1.1 


BD+65 1043 


4615 


5250 


0.0 


<2 


<0.4 


<1.2 


BD+65 1044 


6688 


6500 


0.0 


<3 


<1.6 


<1.5 


BD+74 0718 




5500 


-0.1 


<3 




<1.2 


BD+74 0719 




5500 


-0.1 


<3 




<1.2 


HD 167215 




6500 


-0.2 


43+3 




2.8 


HD 167216 




6500 


-0.2 


34+3 




2.7 
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Table 3 — Continued 



HD/BD 


T eff (J - K) 


T eff (fit) 


[m/H] 


EW (Li I) 


logN(Li)coo 


logN(Li) fit 


BD+69 0993 A 




6500 


0.0 


30±3 




2.7 


BD+69 0993 B 




6500 


0.0 


34±3 




2.7 


BD+08 4386 A 




6000 


0.0 


<2 




<1.6 


BD+08 4386 B 




5500 


0.0 


<3 




<1.2 


BD+03 4428 A 




5500 


-0.2 


<3 




<1.2 


BD+03 4428 B 




5500 


-0.2 


<3 




<1.4 


BD+76 0835 A 




6000 


-0.2 


<2 




<1.5 


BD+76 0835 B 




5750 


-0.2 


<2 




<1.3 


BD+01 4575 A 




5750 


-0.2 


<3 




<1.2 


BD+01 4575 B 




5750 


-0.2 


<3 




<1.6 


BD+11 5033 A 


5162 


5250 


-0.1 


<3 


<1.0 


<1.2 


BD+11 5033 B 


5779 


5000 


-0.1 


<3 


<1.4 


<1.1 


HD 224984 




6000 


-0.1 


35±3 




2.3 


HD 224994 




6000 


-0.1 


33±3 




2.3 



